Charge transfer between epitaxial graphene and silicon carbide 
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We analyse doping of graphene grown on SiC in two models which differ by the source of charge 
transfered to graphene, namely, from SiC surface and from bulk donors. For each of the two 
models, we find the maximum electron density induced in monolayer and bilayer graphene, which 
is determined by the difference between the work function for electrons in pristine graphene and 
donor states on/in SiC, and analyse the responsivity of graphene to the density variation by means 
of electrostatic gates. 



Graphene [T] - monolayer or bilayer of carbon atoms 
with a honeycomb lattice - is a gapless semiconductor, 
which can be used as a current-carrying element in field- 
effect transistors. It has been found that the epitaxial 
graphene grown onto cm-size wafers of the Si-terminated 
face of SiC [5H1] maintains structural integrity over large 
area and demonstrates a relatively high mobility of car- 
riers [SHTj- This makes graphene synthesized on SiC 
(SiC/G) a promising platform to build integrated elec- 
tronic circuits, assuming one can control the carrier den- 
sity in it. For transistor applications, bilayer graphene in 
SiC/G is a particularly interesting material, since inter- 
layer asymmetry (e.g., induced by a transverse electric 
field) opens a minigap in its spectrum [HI E] . 

In this Letter we present a theory of charge transfer 
from SiC to the epitaxial monolayer (MLG) and bilayer 
(BLG) graphene grown on its surface. It has been noticed 
that epitaxial graphene is always substantially n-doped, 
so that use of SiC/G in transistors requires reduction 
of carrier density using gates [3]. The initial doping of 
graphene comes from a combination of bulk donors in SiC 
with the volume density p and surface donor states with 
the sheet density of surface states, 7 (DoS). The charge 
transfer to graphene in a top-gated field-effect transistor 
can be found from solving two coupled equations: 

j[A — 47re^(i(n + Ug) — £_F(n)] + pi = n + Ug, (1) 
A^eF{n) + U + 4TTe'^d{n + ng). (2) 

Equation (nj) accounts for the charge balance, with Ug = 
CVg/e (e > 0) being the areal density of electrons trans- 
ferred to the gate. Here, A j A is the difference between 
the work function of graphene and the work function of 
electrons in the bulk/surface donors in SiC, and ep is the 
Fermi energy in doped graphene. Equation ([2]) describes 
the equilibrium between electrons in graphene and bulk 
donors, with / standing for the depletion layer width in 
SiC, C/ = "^pP being the height of the Shottky barrier 
(x is dielectric constant of SiC), and d - the distance be- 
tween the SiC surface and the middle of graphene layer. 
In the following, we calculate the density n of elec- 
trons (for both MLG and BLG) in two limits: graphene 




FIG. 1: Sketch of a SiC/G-based field effect transistor. Insets 
illustrate relevant part of the MLG and BLG band structure. 



doping dominated by the charge transfer from (a) sur- 
face donors, which corresponds to solving Eq. (fTl) with 
pi — > 0, and (b) bulk donors (e.g. nitrogen), which 
corresponds to solving Eqs. (1J2| with 7 = 0. Charge 
transfer in a more generic situation, with arbitrary p 
and 7, can be assessed by taking the largest of the two 
estimates. Then, we determine the responsivity factor, 
r = —dn/dug, which characterizes the ability to control 
the carrier density using external gates: r — > 1 would be 
optimal for transistor operation of gated SiC/G devices 
whereas r <^ 1 would indicate that transistor operation 
is impossible. 



Monolayer graphene (MLG) has linear spectrum 
£±(p) — i'^Pj in the two valleys, corresponding to the 
non-equivalent corners K and K' of the hexagonal Bril- 
louin zone, so that epin) = fivy/'r:n (we take into account 
both valley and spin degeneracy of the electron states) . 



In the limit (a) we find that the carrier density is 
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(3) 



The initial density of electrons in graphene is described 
by Eq. (|3]) (n^ = 0) with two characteristic regimes, 



A7, 7 < 7f ; 



(4) 



where n^ and 7^ are, respectively, the saturation value 
for the carrier density in SiC/G and the crossover value of 
DoS of donors on the SiC surface at which n(7) saturates: 



lA 



A^A 



{'^^ 



lA 
Id 



(l + y^ 



lA 
Id 



(5) 



The responsivity of MLG in SiC/G is, then, 
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(6) 



In the limit (b) we find n{p) by solving Eqs. ( 1|2 ) 
numerically. The numerical solution shown in Fig. 121 
(d) interpolates between the regimes of weak and strong 
graphene doping: 
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The responsivity, r|„ =0 of MLG in SiC/G in the limit (b) 
can be described using Eq. (p|, but with ja replaced by 
jA and the upper/lower limits corresponding to p <C pi 
and p^ pi, respectively. 

Bilayer graphene (BLG) has the spectrum [Sj with 
conduction (+) and valence (-) bands, which, in the vicin- 
ity of the Brillouin zone corners, are described by 



where, ei is the interlayer coupling and s distinguishes 
between pairs of degenerate (s = —1) and split (s = 




FIG. 2: Comparison between charge transfer from SiC to 
MLG and BLG, with A measured in units of split-band energy 
in BLG, ei = 0.4 eV, d = 0.3 nm for MLG and d = 0.5 nm 
for BLG. (a) Electron concentration in graphene n dominated 
by charge transfer from surface states, (b) Values 7* of the 
surface DoS at which n{'y) saturates, (c) Saturation density 
value as a function of n«, in units of ni. (d) Electron bulk 
donors density p for 7 = 0. 



1) bands. The Fermi level in n-doped BLG with sheet 
electron density n is determined as 



In the limit (b), when n is determined by charge trans- 
fer from bulk donors in SiC, 
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2 , n <ni, 
n > ni, 



where m = 2£f /(ttS-^u^) w 3 • 10" cm"^. 
In the limit (a), this gives 
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for n (7) < ni and Ug — 0. Here, 71 = ei/Trh v and 
71 = 71 + 7a/2. For larger densities, n^' > ni, 
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which resembles Eq. (|3| for MLG, but with 7^ — >■ 27^1. 

Similarly to MLG, the density in BLG on SiC saturates 
upon the increase of surface DoS of donors. The crossover 
to the saturated density, nl^ {A) — A'yl^{A), occurs at 
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The dependence of n^^{'-f), j^^ {A) and nl^ (A) on the 
relative size of the band splitting ei, and the graphene- 
surface donors work function A is shown in Fig. [2] (a-c). 
Responsivity of the BLG to the gate voltage is high or 
low, depending on whether the saturation regime for the 
carrier density is reached, or not. For 7^7, 
For 7 > 7^ 
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(A), and responsivity r w 1 of BLG re- 

responsivity 
ISI), with 7A 



r.11 



with h 

quires that p <C pi^ , whereas for p ^ pi 
is described by the same limits as in Eq. 
replaced by 7^. 

Independently of the number of layers, the gate volt- 
age V* needed to reach the neutrality point in graphene 
controlled by the top gate with mutial capacitance C is 



V* w - max 
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In summary, we calculated charge transfer from SiC to 
epitaxial graphene. In the case when the charge trans- 
fer is dominated by donors on the surface of SiC with 
A - leV or donors in the bulk of SiC with A - leV, [H] . 
we estimate that the saturation density of n-type doping 
of MLG is lO^'^cm"^, which corresponds to ep ~ 0.4eV 
(for d « 0.3nm). This value of carrier density occurs 
when the donor volume density is p > p^ ~ 10^^ cm~^ 
(we use X ~ 10 for 6H SiC [13]) or the surface states have 
DoS 7 > 7^^ ^ 1 • 10^'^ cm~^eV~^. For lesser doping of 
SiC, 7 < 7* and p < p^, o ne should use the larger of the 
estimates from Eqs. (4J7) This can be compared to the 
data reported in the recent studies of epitaxial graphene 
indicating a substantial initial level of n-type doping of 
SiC/G, very often [4] as high as lO^^cm"^. However, 
some particular growth processes produce SiC/G with a 
much lower doping level [51 [TT] , indicating that efficient 
annealing of donors on and near the SiC surface is pos- 
sible. 
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